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Abstract. - We describe the Fermionic Casimir effect in the case of two spherical superfluid 
C^ ■ scatterers immersed in a normal Fermi system. It is shown that due to the focusing property 

of Andreev reflection this new Casimir-like energy is significantly enhanced when compared to 
the case of normal scatterers with specular reflection. 

o 

Introduction. - If two or more voids (or hard-wall scatterers) are formed in a normal 

CN| , uniform Fermi system, effective Casimir-like interactions between such voids appear [1,2]. For 

^ ' example, two spherical voids of radius a with distance between their centers _R in a Fermi 

ir-s ' medium at zero temperature and with Fermi wave vector fcp would interact with each other 

P^ ' with the effective potential energy 

O [ EaW . -.,JJ^^-^^^M2,,(R - 2a)l , (1) 

where t'dog is the spin-isospin degeneracy factor. Here ji{x) is the spherical Bessel function 
Ci I and m the fermion mass. There is an illuminating explanation and easy derivation of this 
expression based on semi-classical physics — periodic orbits theory in particular. One can 
show that the energy density of an infinite and otherwise homogeneous Fermi system is affected 
by the presence of the two spherical voids and that the main contribution to this energy density 
Q I change is determined by the "almost standing wave" created between the two scatterers. As 
O ' was shown in refs. [1,2], this argument can be extended, in principle, to an arbitrary number 
of scatterers of arbitrary shapes. Volovik has arrived at similar conclusions by studying 1- 
dimensional systems [3] and he refers to this energy as the Mesoscopic Casimir energy or 
force. In his discussion of this new type of Casimir energy, Volovik stresses that this energy 
3 ■ is governed by the "microscopic trans-Planckian physics," unlike the usual Casimir energy, 
which is determined by the long-wave physics. 

The Fermionic Casimir energy appears because of the Dirichlet boundary conditions at 
the voids, which semi-classically correspond to specular reflections. If, however, the boundary 
separates a normal Fermi system from a superfluid Fermi system, the specular reflection is 
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Fig. 1 - Two superfluid grains Si and S2 immersed in a normal Fermi gas N, is the angle at 
center specifying the position of the emerging ray on the first disk/sphere and 9 the angle made by 
its momentum with the normal. 



replaced with Andreev reflection [4]. A particle at the Fermi energy in the normal part of 
the Fermi system, impinging on the boundary between the normal and superfluid system, 
is retro-reflected as a hole, irrespective of the shape of the boundary. This rather simple 
but amazing fact has apparently quite unexpected consequences. Let us imagine now that 
we have two spherical superfluid grains embedded in a normal Fermi system. One might 
then naively expect that there is an infinite number of stable (classical) periodic trajectories 
between two such objects. This situation has to be contrasted with the case of specular 
reflection, when only one unstable trajectory between the "noses" of the two spheres exist. 
The latter trajectory is responsible for the change in the density of states, which ultimately 
results in the Fermionic Casimir energy discussed in refs. [1-3]. One can now expect that an 
infinite number of stable periodic trajectories would lead to a very large change in the density 
of states if two superfluid grains are embedded into a normal Fermi system. From this point 
of view the fact that two spherical superfluid scatterers do not generate a chaotic dynamics 
is initially quite a surprise. Our goal here is to determine the character and magnitude of 
the Fermionic Casimir energy in the case of two spherical superfluid scatterers embedded in a 
normal Fermi system. The result we obtain can easily be generalized to scatterers of arbitrary 
shapes and to any number of scatterers in a leading order approximation, which we shall 
define more carefully below. We would like to mention here that even though other aspects 
of the physics of superfiuid grains of various sizes embedded in a normal Fermi system have 
been addressed before (see refs. [5,6] and references therein), the existence and nature of the 
Fermionic Casimir energy was never considered. Various other aspects of the quantum chaos 
in so-called Andreev billiards/stadiums have been considered by many authors [7] , but spatial 
configurations of the type we shall study here have not. Our interest stems partially from 
the physics of neutron stars, where such structures are expected to exist [1,8]. Moreover, the 
possible relevance of Andreev reflection was mentioned in connection with the appearance of 
QCD-color-superfluid "grains" in the deconflned quark-gluon phase in the core of neutron [9] 
or of quark stars, and also in connection with the observed cosmological asymmetry between 
matter and anti-matter and the origin of dark matter [10]. 

Two superfluid grains. - To set the stage, let us consider the geometry illustrated in fig.Q 
In the Andreev approximation any straight line joining two arbitrary points on the surface of 
the two superfluid grains leads to a classical periodic orbit. A particle (p) leaving the surface 
of one grain, hits the second grain somewhere and it is reflected exactly backward as a hole 
(h) and it thus returns to the initial point on the flrst grain, where it now transforms, after 
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Fig. 2 - The Poincare surface of section for the case of two supcrfluid circular/spherical grains 
embedded in a normal Fermi gas (in two dimensions) , separated by a distance four times larger then 
their respective radii {R=Qa). The ratio of the particle and hole momenta is here fcp/A:h = 1.5. Only 
stable/elliptic trajectories parameterized by the angle (j) (in degrees), where the particle originates on 
the surface of the disk Si, and by the angle (j> — 6 (in degrees) between the momentum of the particle 
and a parallel of the center-to-center axis are shown. The relatively small values oi (j) — 6 show that 
most trajectories run in a tube around the line joining the two centers. 

another Andreev reflection, into a particle. This is strictly correct if the particle has exactly 
the Fermi energy. Then the energy of both particle and hole states are equal. Otherwise 
the relation between the angles of the incoming particle/hole and the reflected hole/particle 
resembles Snell's law: 



2 2 



/Cpsin0p==fchsin0h, ^ + ^ = 2Ai = fc|, (2) 



where /x is the Fermi energy and kp > kp , kh < kp denote the fermion momentum in a particle 
state and a hole state, respectively (we shall use units: h^ = m = 1). Note that the energies 
of the particle and hole states are symmetric with respect to the Fermi energy. Moreover, 
the incident and reflected rays are on the same side of the normal to the surface at the point 
of incidence, as opposed to the case of specular reflection. Consequently, one can show that, 
unlike in the case of normal reflection, the disks/spheres now have a focusing effect. 

Instead of trying to seek various possible periodic trajectories we simply generated a 
Poincare surface of section. Only those orbits are stable for which the separation between 
the centers of the grains satisfies the following constraint, derived by requiring that the eigen- 
values of the monodromy matrix have to be inside or on the unit circle in the complex plane: 

— < — ^ , (3) 

where we have assumed that \kp — kh\/kp ^ 1. If this constraint is not satisfied the orbits have 
a hyperbolic/unstable character. In the case of equality the eigenvalues of the monodromy 
matrix are equal to (minus) one, the orbits are parabolic/marginally stable in character and 
elliptic/stable otherwise. 

Exact expression of the density of states. - The structure of the classical phase space 
briefly described here would apparently imply a rather complex quantum mechanical dynam- 
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ics. As we shall show, however, the most important features of the density of states in such 
systems can be understood and described quite accurately in terms of a simplest periodic 
orbit and its repetitions. In order to determine the change in the local density of states we 
shall construct explicitly the scattering states and determine the phase shifts by extending 
the formalism used in case of normal systems [2, 11, 12] to Bogoliubov-de Gennes equations. 



A*(r) 



A(r) 
-(T-/i) 



u{r) 
v(r) 



E 



u{r) 
v{r) 



(4) 



where T is the kinetic energy operator, E is the quasi-particle energy, u{r) and v{r) are 
the particle and hole spatial components of the quasi-particle wave function and A(r) is the 
pairing field, assumed to be constant inside the two scatterers, i.e. |A(r)| = A, and vanishing 
outside. For the sake of simplicity of the analysis we assume that the mean field vanishes 
everywhere and that the two scatterers have the same radius a. Physically it makes sense 
to discuss many fermion systems for which fcpa 3> 1. One can distinguish further between 
small and large superfluid grains, depending on whether the coherence length ^ — kp/nA is 
smaller or larger than the radius of the grain a. The pairing fields in the two grains need not 
to have the same phase, which for the sake of simplicity of the argument we consider to be 
constant over the grain. As in ref. [2], we shall use Krein's formula [13], generalized to the 
case of Bogoliubov-de Genes equations. Krein's formula expresses the change in the density 
of states through the phase shifts, and represents a generalization of a much simpler formula 
due to Beth and Uhlenbeck [14]. We thus obtain in the 2-dimensional case (the 3-dimensional 
case has a very similar structure): 



6g{E) 



dSN{E) 
dE 



dN{E) 



-lmlndet(M) , 
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(5) 
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r(i) 



r(2) 



r(i) 



r(2) 



i^p(a, E, fi, A) H];^>_^,{k^R) f™,(a, E, m, A)e^*-i/„_„ 



,{KR) 



(2) 



t^^{a, i?,M, A)e-*-H,^„^i„,(fcpi?) t™ («, E, f,, A) <i„,(fchi?) 



with kl = 2{fi 



E), k^ = 2(^ — E) and where to, to-' = 0, ±1, ±2, • • • are angular momenta. 
The angle 0a is the relative phase of the pairing fields in the two grains and II is the identity 
matrix. The matrix elements t™ = tVl{a,E,iJ,,A), . . . are equal to half of the T-matrix 
elements for the scattering of a particle/hole of angular momentum m off a single spherical 
superfluid grain (centered at the origin of the coordinate system), and the ordinary Hankel 
functions H^_^,{kpR), . . . determine the propagation between the grains. 

Asymptotic approximations. - Even though these formulas are already suitable to per- 
form explicit calculations, one can take advantage of the fact that k-pR ^ 1 and obtain 
signiflcantly simpler expressions by retaining only the leading term in an asymptotic expan- 
sion of the corresponding Hankel functions. One thus arrives at the following much more 
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Fig. 3 - The structure functions Tpp, Thh and Tph for the case of two identical disks with /i — 200 /a^ 
and A = 50/a^ and (^a = 0. The sohd hues are the real parts of the corresponding quantities and 
the dashed lines the imaginary parts, respectively. Notice the difference in the magnitudes. 



transparent expressions 

2 



^ph = (X^^pii) cos(/)A = |rph|exp(i0ph)cos0A, 

^pp = (E Cp(-1)'")' ^ I^PPI e^P(*'^Pp) ' ^hh = (E Chl-l)")' ^ iThhl exp(z0hh) 



^^^^^ « &^i^sin[(fcp-A:h)i? + 0ph]-^^cos(2fcpi?- 



TT^^/Ch/Cpi? 

, 2|Thh' 



71-2 fcpi? 



-^pp; 



7r2A:hi? 



cos(2A:hi? - (/)hh) • 



(6) 



The structure functions Tp^, Tpp and T^h characterize the properties of the scatterers and 
are independent of the separation R between them. In fig. |21 we show some typical behavior 
of these functions with the quasi-particle energy E. Note the great difference in size between 
the normal (pp, hh) and Andreev (ph) structure functions, especially when E < A. 

Even though these structure functions are relatively easy to evaluate numerically in the 
case of almost any conceivable size grains, it is extremely instructive to go one step further 
and derive an approximate expression for these structure functions themselves. An analysis of 
the structure functions (omitted here) eventually leads us to the following significantly simpler 
expression for the change in the density of states valid for i5 ^ |A|: 



SN{E) = -Im ■ 



a niin{kp ^k]^) 



du exp 



a min(A:p,A;h) 



i{kp - fch)t^^ 



„i(fcp-fch)(fl-2a) 
TT^i^^fcpfch 



cos((/)a)- (7) 



One can easily recognize the argument of the second exponential as the classical action of the 
shortest particle-hole periodic orbit. The hole momentum enters with a minus sign since it 
is anti-parallel to the hole velocity. The squared integral would become a product of similar 
factors in the case of two grains of different radii. For £'=0 this overall factor becomes 
Ami 



,2 
"max 



Ak^a^ i.e. the squared number of partial waves undergoing Andreev reflection. 
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Fig. 4 - The density of states variation SN{E) evaluated exactly (solid line), in the approximation 
of eq. ||SJ (dashed line) and under a further asymptotic expansion of the cylindrical functions in the 
structure functions (dotted line) for two identical disks with R—Qa, jj. — 200/a^, A=50/a^ and (f>A~0. 



Results and conclusions. ~ In fig. ^ we compare the change in the density of states 
evaluated exactly (see eq. (O) vifith the simplified expression of eq. © and with a further 
asymptotic approximation of the structure functions on which eq. (Q is based. It is now 
obvious that the oscillatory character of the change in the density of states is fully accounted 
for by the shortest periodic orbit, quantized as expected by 



{kp - A;h)(i? - 2a) - 2 arccos (E/A) = iirn , n = 0, ±1, ±2, 



(8) 



The maxima in Sg{E) = dSN{E)/dE correspond exactly to those energies for which the above 
condition is fulfilled. At each Andreev refiection the reflected wave acquires an additional 
phase arccos(ii^/A), leading to the overall minus sign in eq. 0, valid for the case E ^ A. 

By using the exact expression Q or the simplifications ((HJ or Q, one can now easily 
evaluate the interaction energy of spherical superfluid grains in the Fermi sea. 



/■OO 

Esc{R)^i^dce dEE 

Jo 



dSN{E) 
d^ 



— — l^dcg 



dESN{E). 



(9) 



as function of their separation R. This Casimir energy is completely dominated by the states 
with _E < A. We thus obtain, by numerical integration, for two spherical three-dimensional 
grains 



£^SC (R) ~ i'dcg 



h^k^. 



kpa 



2m 27ri?2(i?-2a) 



X . 



(10) 



where X = 0(1) and positive and very weakly dependent on separation. By comparing with 
eq. (^, one can show that Esc{R) ^ l^c(-R)| for the typical case kpa ^ 1. This is easy to 
accept now: in the case of two superfluid grains, embedded in a normal Fermi system, there are 
many "bound states" (see fig. |2) due to the appearance of a large number of periodic orbits, 
whereas there exists only a single unstable periodic orbit between two hard spheres [1,2]. 

In summary, we have constructed both exact and simplified expressions for the calculation 
of the density of states and the Fermionic Casimir energy in the novel case of two spherical 
superfluid scatterers immersed in a normal Fermi system. We have furthermore shown and 
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explained that these quantities are significantly enhanced in comparison to the normal case 
where the obstacles are scatterers with hard-wall reflection conditions. Since the change in the 
density of states is dominated by the Andreev reflection, which is almost exactly backward, 
the case of many grains reduces simply to the sum over all unobstructed pairs of grains. In 
particular, one can expect that the energetics in portions of the neutron stars would be greatly 
affected when the Fermionic Casimir energy is taken into account — more so than was the 
case when only normal systems were considered [1,8]. Generalizations of the present results 
to more complicated geometries appear to be straightforward. 
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